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Abstract 

We consider current precision electroweak data, Z' searches and dark matter constraints 
and analyse their implications for an extension of the SM that includes an extra U(l)' 
massive gauge boson and a particular hidden sector ("hidden valley") with a confining 
(QCD-like) gauge group. The constraints on the Z' with arbitrary Z — Z' kinetic mixing 
coming from direct searches and precision tests of the Standard Model are analysed and 
shown to lead to a lower limit of 800 GeV on its mass. Renormalisable interactions 
involving the Z' probe the physics of the hidden valley sector which contains a pseudoscalar 
dark matter candidate. We find that dark matter constraints place an upper bound on 
the mass of the Z' of O(10) TeV. A TeV mass scale is needed for the hidden valley states, 
and the Sommerfeld factor for p-wave dark matter annihilation is found significantly to 
suppress the allowed parameter space of the model. 
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1 Introduction 

The Standard Model (SM) has so far proved remarkably successful. With the Large Hadron 
Collider (LHC) about to start operating, confirmation of its Higgs sector or searches for new 
physics beyond the SM enter a new, interesting stage, when theoretical ideas and models 
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advocated over the past three decades face the test of new experimental data. In this paper 
we consider current precision electroweak data and dark matter constraints and analyse their 
implications for a rather simple extension of the SM. This consists of the SM and a particular 
hidden sector that contains matter which interacts under a new confining gauge group, with 
an extra U(l)' massive boson that couples to SM and hidden sector matter. It was partly 
motivated by novel LHC signatures and search strategies specific to this type of model. Further 
interest in this class of models has also been triggered by attempts to explain the excesses in 
the positron spectrum measured by PAMELA [lj and ATIC [2]. 

Hidden sectors are a common presence in model building beyond the SM, in models with 
branes in extra dimensions and in string and supersymmetric models, with many implications 
such as, for example, supersymmetry breaking [9]. The physics of the hidden sector can in 
principle be probed by the mediating particles communicating between the hidden and visible 
sectors. While gravitational interactions are the usual standard mediators, renormalisable 
operators involving messenger fields such as new gauge bosons, scalars, neutralinos or mes- 
senger fermions may also be present. For this latter case, a special class of hidden sectors 
is that advocated in [3J, also called "hidden valley" (HV), where the hidden sector contains 
some light states in comparison to the energy scale of the mediator. If these light states are 
able to decay to SM states, the dynamics of the HV models can then produce some novel 
phenomenology [3j 01 [5j El [8], possibly relevant to the LHC. To select search strategies 
for such new physics, the restrictions that exist on these models from electroweak and dark 
matter constraints need to be well understood, and this is the main task of this paper. 

In this paper we present the constraints on a "hidden valley" extension of the SM with a 
new neutral gauge boson (Z') dominantly responsible for mediation between the SM and HV 
sectors. In Section[2]we present the hidden valley model considered, based on the construction 
of Strassler (5J, and define the parameters that will be constrained by current experimental 
data. In this model the HV matter is confined by a strong (QCD-like) interaction in the hidden 
sector. The physics of the confining theory is discussed and the spectrum and stability of the 
HV states determined. In Section [3J we review the constraints on the Z' from electroweak 
precision data (EWPD) and direct searches, leading to a lower mass bound for the mediator. 
The discussion and conclusions of this section are more general and independent of the physics 
of the hidden sector. The numerical limits presented allow for the possibility of arbitrarily 
large kinetic mixing between the hypercharge and U(l)' gauge bosons. So far as we are aware, 
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these are new results and apply to the case that the SM matter have 0(1) U(l)' charges. In 
Section 0J the consequences of a dark matter candidate in the hidden sector are considered, 
leading to constraints on the HV dark matter candidate. This will produce an upper mass 
bound on the Z' mediator. We demonstrate that the lower and upper mass bounds can be 
simultaneously satisfied in a finite region of parameter space that we identify in detail. 

The Appendix provides some calculations and results used in the text. In Appendix [X] 
additional details on neutrinos masses are provided, while in Appendix [B] the gauge boson 
eigenstates are found after the addition of the local U (1)' symmetry to the SM. In Appendix ICl 
we list the perturbative gauge boson partial decay rates and spin-0 particle annihilation cross 
sections used in Section HI In Appendix [Dl we discuss non-perturbative corrections to the 
matrix element for slow moving particles, and present the Sommerfeld factor for p-wave mul- 
tiparticle states in the presence of Yukawa interactions used for the dark matter relic density 
calculations. 

2 A "hidden valley" extension of the Standard Model 

Here we construct a simple extension of the SM based on the construction in [3], which 
can lead to an interesting non-standard phenomenology. The choice of hidden valley matter 
content and symmetries is somewhat arbitrary. In addition to the non-supersymmetric SM 
particle content with right-handed neutrinos (qi,Ui, di,^, e i , N i} H), four new Weyl fermions 
(U L R , D L R ) are introduced that make up the hidden sector and a new scalar <j) is added. 
The new states are singlets under the SM gauge group G SM . An extra U(l)' gauge group 
is also assumed under which both SM fields and hidden sector ones are charged, and so by 
construction the associated gauge boson mediates between the SM and HV sectors. 
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Table 1: Charges of the fields under the U{\)' that extends G, 



A simple way to ensure the cancellation of anomalies associated with the U(l)' is to choose 
the U(l)' charges for the SM fields to be those of the family independent U(l) x subgroup 
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of a 5*0(10) GUT, as shown in Tabled! As discussed in Appendix [A] this model may be 
considered a low energy realisation of a theory with 50(10) symmetry. However, we do not 
restrict the analysis to be consistent with 50(10) unification, but this may be recovered as 
a special case. It follows with the particle content and £7(1)' charge assignment in Table [T] 
that the "un-unified" HV model is still gauge anomaly free. Other U(l)' extensions are 
possible and will bring different numerical limits on the model parameters compared to those 
presented in sections [3] and |H However, as we shall see, the qualitative features of the allowed 
parameter space found are universal. The vacuum expectation value (vev) of 4> is responsible 
for generating a mass for the Z' which we will find is of 0(TeV). To allow for right handed 
neutrinos masses we choose Q', = 2. 

One constraint on the hidden valley model which must be satisfied is the upper bound on 
the number of relativistic degrees of freedom in thermal equilibrium at the time of Big Bang 
Nucleosynthesis (BBN). This can be satisfied by ensuring that there is a sufficient mass gap for 
the new states that this model introduces. In the hidden sector, interactions other than £7(1)' 
that are restricted to this sector can also exist. We assume as in [3] that there is a confining 
SU(N V ) QCD-like interaction with the HV fermions in the fundamental representation and 
with a confinement scale A HV . If A HV > A QCD , a sufficient mass gap would be provided to satisfy 
the BBN constraint. This hierarchy can naturally set the mass scale of the HV states to be 
significantly different from that of SM states, using arguments similar to those in technicolour 
models. The spectrum of bound states is discussed in Section [231 

To allow current masses for the HV fermions via the Higgs mechanism, the condition 
q + + <7_ = — Q', must be satisfied. The field eft can then also act as a mediator between 
the SM and HV sector, via interactions with the SM Higgs of the form \H\ 2 \4>\ 2 . This 
term produces mass mixing terms between the scalars as well, after breaking the £7(1)' and 
electroweak symmetry. The use of a scalar mediator (also known as the "Higgs portal") has 
received recent attention in [TU [T21 [T31 EH EH]. In this paper we expect that mediation via 
the Higgs portal is negligible due to small Yukawa couplings, and the large mass for the <f>, 
with the resultant small Higgs-0 mixing, making the Z' interactions dominant. Indeed the 
Higgs portal can be ignored completely in this model if the HV fermions current masses are 
entirely generated from instanton effects after spontaneous chiral symmetry breaking [15J. 
Alternatively, if q+ + g_ = 0, mass terms for the HV fermions are unprotected by gauge 
symmetries and they do not couple directly to the fundamental scalar fields. 
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2.1 U(l)' physics. 

The addition of the U(l)' gauge boson to the SM introduces interesting physics by itself 
[1Q\ 117]. The discussion here, and in Section [3] for the limits on the Z' from corrections to SM 
processes, is independent of the details of the hidden sector physics. In the most general case, 
the Lagrangian can contain a term that leads to kinetic mixing between the Abelian gauge 
fields U(1) Y and 17(1)': 



C 



1 



gauge kinetic 



[Fy F Y ^ + F'^ F'^ + 2 sin x F^Fff) + non-Abelian kinetic terms. (1) 



The presence of kinetic mixing and mass mixing terms changes the currents coupling to 
the gauge fields from that expected in the SM. Since the new scalar, <fi, is a singlet under 
G SM , the only unbroken Abelian symmetry after f7(l)' and electroweak symmetry breaking 
(EWSB) is that associated with electric charge conservation, generating the usual structure 
for electromagnetism. Thus, the physical photon has the same couplings as in the SM, and 
there is no milli-(electrically)-charged matter in this model at the tree level. After EWSB, an 
F'F^m term is generated by loops. For U(l)' couplings in the perturbative regime, this term 
is negligible for the following analysis. 

The currents Jj^zz'i which include the relevant coupling constants, associated with the 
canonically normalised gauge boson mass and kinetic eigenstates of the model are related to 
SM and £7(1)' currents </|m Z , J' as follows (see appendix A for details): 
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(3) 



(4) 



The £ angle accounts for the rotation necessary to diagonalise the mass matrix, after kinetic 
mixing has been removed by a field redefinition. Note that if £ = 0, the coupling to the physical 
Z is identical to the SM case, but the Z' coupling still varies with the kinetic mixing x- It is 
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therefore not necessary for any SM particles to be charged under the U(l)' for the physical 

Z 1 to mediate between the SM and HV sectors as this could be achieved by the presence of 

kinetic mixing alone. In eq ([2]), the Weinberg angle is defined as W = arctan(g y / g w ) , and 

Q' H in eq ^ is the U(l)' charge of the SM Higgs field. The tree-level mass of the physical Z 

deviates from the SM prediction: 

2 2 rr&r sec 2£ sec 2 W + mi(l — sec 20 

m 7 = (5) 

z l + sec2C K ' 

The presence of a local U(l)' symmetry thus affects "low energy" electroweak physics whose 
parameters have been measured very accurately. The impact of experimental measurements 
on the mixing parameters, Z' mass, U(l)' coupling constant and matter charges are analysed 
in detail in Section El 

2.2 Physics of a QCD-like interaction in the hidden valley 

Let us now address the details of the physics of the hidden sector. We consider the case 

that the HV matter is confined by the strong dynamics of an unbroken gauge symmetry in 

the hidden sector. We review the behaviour of such QCD-like theories where the elementary 

matter fields are vector-like pairs of fermions and note some possible differences with SM 

QCD phenomenology. Here we focus on the lightest composite states since, if stable, these 

could be the dominant dark matter candidates. The experimental constraints on the hidden 

valley model with a dark matter candidate are discussed in Section HI 

The Lagrangian with elementary fermions in the fundamental representation of a confining 

SU(N V ) symmetry (N v > 2) includes the terms: 

1 q 2 ' 
4 y ^ y 32 tt 2 



p _ pC pCflU y P HV pC pC^lV ff>\ 

^HV gauge A » ' fll/^ 1 on 2 ^ ilV^ V° ) 



£ H vf Crm io„ = 9 L 0* L + V R il}* R -* R M* L -* L Mt* R (7) 

where Q is the field strength of the HV-confined gauge field, and chiral multiplets {^>) have 
been formed out of the left- and right-handed Weyl fermions (U,D). Axial rotations can be 
used to redefine the valley-quark (or "v-quark") fields, which allows phases to be removed from 
the v-quark mass matrix. However, a U(1) A transformation generates non-zero surface terms 
which contribute to the QQ term [20]. The effect of choosing a basis where the transformed 
mass matrix, M', is real and positive semi-definite, is to shift the # HV parameter: 

0hv -»• H v + argdetM = nv (8) 
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While in the SM sector, |#q CD | < 10~ 9 , in the present case nothing prevents # HV from taking 
values of order O(l). This term explicitly violates CP symmetry unless sin 9 = 0. 

The dynamics of the confined matter can be modelled using an effective theory built 
to possess the conserved and approximate symmetries, as detailed in the following. The 
kinetic terms of £ HV fermion are invariant under independent transformations of the left- and 
right-handed multiplets. Thus, there is a global U(N) L x U(N) R symmetry, where N is the 
number of HV fermion flavours. The mass terms explicitly break the chiral symmetries down 
to a vector subgroup. If all flavours have a common current mass, the final symmetry is 
U{N) V . If the current masses are not degenerate, the final symmetry is [U(1) V ] N . 

If there were no fermion masses, the axial symmetry would be spontaneously broken by 
the v-quark condensate whilst still preserving the vector symmetries. Thus, if there are 
N' fermions whose current masses are much smaller than the confinement scale, then it 
is a good approximation to assume that the chiral symmetry of the kinetic term is first 
explicitly broken to U(N') L x U(N') R x [U(1) V ]^ N ~ N '\ and then spontaneously broken to 
U{N') V x [U(1) V ]^ N ~ N \ The strong interaction is flavour blind in the limit of zero current 
masses, so does not spontaneously break the vector symmetry further. The pseudo-Goldstone 
bosons associated with the spontaneous breaking of the approximate axial symmetry are the 
analogues of the pions, kaons, eta and eta primed mesons in QCD. 

It is convenient to use chiral perturbation theory which is an expansion in terms of the 
explicit symmetry breaking parameter, M/A HV , to model the pseudo-Goldstone boson com- 
posite states. The degrees of freedom are parametrised by a field matrix, E, which is a function 
of the ./V' 2 pseudo-Goldstone bosons: 

= ( 9 ) 

where v accounts for the magnitude of the quark condensate and has mass dimension one, 
and E is an N' x N' matrix. Under a U(N') L x U(N') R transformation, the fields transform 
as follows: 

LEi? f (10) 



n = 

K = 
E' = 
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An exponential representation for E is chosen here, defined using the decay constant, f n : 



2 * n ^, U^-^=l+n a T a . (11) 



The T a generators are normalised so that Tr[T a T b ] = S ab /2. The pseudo-Goldstone bosons 
found after spontaneous axial symmetry breaking are pseudoscalars. The effective Lagrangian 
for E, to leading order in M and momentum, is given by [18J : 



f 2 

L ~ 4 



Tr 



D^D^] + (2^Tr [EM] + h.c) + ^ ( ^ - 1 ] Tr [E^E] 



(12) 



where /i = 2v 3 / f%. Although the axial C/(l) symmetry (L = i?t) of the kinetic terms is 
classically conserved in the limit M — > 0, it is always broken by the triangle anomaly at the 
loop level. The effect of the anomaly in a SU(N V ) confining theory is vanishing in the limit 
N v — > oo. To describe physics at finite N v , perturbations can be added with increasing powers 
of 1/N V |19j . This is important for correctly describing the r/ pseudo-Goldstone boson. 

The triangle anomaly also controls the QQ term. It has been shown by 't Hooft [20j that 
this term can be rewritten as a determinantal interaction of the quark fields, det^g]. The 
leading order correction in 1/N V to the effective Lagrangian, which explicitly breaks the U{1) A 
symmetry, is then given below: 

5 a C oc e ie Hvdet[E] + h.c. (13) 

The vacuum alignment of the quark condensate in the flavour space, (E), is selected where the 
potential of the effective theory is minimised. Once this is found, the fields can be expanded 
about their vevs to determine physical parameters such as the masses and couplings. 

Before proceeding further, along this line, let us consider the nature of the dark matter 
candidates. As discussed later in Section 14.21 the flavoured "valley-pions" (or "v-pions," 
-n"^) are the dark matter candidates, which is possible since a flavour symmetry can protect 
them from decay. The pion analogues are the lightest states in this effective theory, with 
the 7r^ (= 7r 3 ) being flavour neutral, and 7r^ (oc tt 1 =p Z7r 2 ) flavoured. The ±,0 indices of the 
v-pions refer to their flavour isospin as they are electrically neutral. It will be shown that the 
dominant annihilation mechanism of the stable v-pions is via the Z' when A HV > 30 m n or 
if there is a fractional mass splitting of 0(10%) between the (heavier) 7r^ and the 7r^. In the 
remainder of this section we show that such a mass splitting is indeed possible. 
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To first order in M, the v-pion masses are degenerate (neglecting any mass mixing effects): 



m 



2 

vp 



Li W m 2 a + m\ + 2 m a m b cos 6> HV (14) 



where m a 6 are the magnitudes of the HV fermion current masses. This degeneracy is lifted 
by 0(M 2 ) termsS 



■J IT 



2 r 



= -f 



( Cl Tr[SM] 2 + c 3 Tr[SMSM] + h.c.) + 2c 2 |Tr[EM]| 2 . (15) 



To find the masses of the v-pions, only the two lightest flavours need to be considered when 
mass mixing can be neglected. For an SU(2) flavour symmetry, one of the terms in eq (|15p is 
redundant, since it can be transformed into one of the other terms using a determinant that 
is not dynamical: 

Tr[A 2 } = Tr[A] 2 - 2Det[A] (A is 2 x 2 matrix) (16) 

Using this we may absorb the c 3 term into c 1 . With these considerations, one obtains the 
v-pion masses to 0(M 2 ) below: 

2 _ 2 / . s m t + m t + 4 m a"ife (ml + m b) cos e " HV + ml ml [2 + 4 cos(2^ HV )] 

m ± — m + (Cj + c 2 j — = 

Wv m£ + + 2 m a m b cos # HV 

<c « ^ + J^Mi^L (it) 

u m£ + + 2m a m b cos fl HV 



which are correct up to O 



CiSm hv v^^v j 

contributions from instanton effects to pseudo-Goldstone boson masses at 0(M N '- X ), see [21] 
and references therein. This however does not affect the mass splitting of the v-pions. 

The best fit values for the coefficients of the higher order terms in QCD [22] determines 
c qcd = 10 _ 3 ^2/j2)qcd x [_o.6 ±0.3, 0.2 ±0.3, (0.9 ± 0.3])(^/3)] i , (N v = 3 for QCD), where 
the N v scaling for follows in the large N v limit [18J. In our case, if the hidden sector physics 
can be approximated by a scaled up version of QCD, using eq (I14p and (I17p with (9 HV ~ 0, a 
10% mass splitting for the valley-pions would be found if (m ab /Li) HV > O(0. 1). This can be 
satisfied with a suitable choice of m ab (< A HV ). If 9 HV ~ tt, the 7r^ fields become massless 
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terms. For a discussion on the 



x The degeneracy is also lifted if the bosons interact non-universally with other fields. 
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in the limit of degenerate quark current masses. However, the 7r„ field remains massive, so it 
can be very easy to find a 10% mass splitting in this regime. As mentioned above, this mass 
splitting plays an important role in the discussion of the dark matter experimental constraints 
on the model, and we return to these in Section SJ 



3 Collider constraints on the Z' boson 

In the following we examine the experimental constraints on the parameters describing the 
U(l)' extension of the SM presented in Section [21 and investigate in particular the bounds on 
the mass of Z'. To this purpose we consider precision measurements at colliders and direct 
searches for new neutral gauge bosons. For a Z' mass of 0(100) GeV or less, additional 
constraints are relevant from determination of the muon anomalous magnetic moment [23j. 
However, for the model considered this does not further restrict the limits presented. 

3.1 Electroweak precision data 

In the following we use a model independent formalism [24J, developed to quantify the devi- 
ation of physics from the Standard Model predictions. The parameters p, x and y introduced 
below define the low energy effective Lagrangian with the physical and Z integrated out: 

AG 

£ c« = F 9Q [sec 2 e w J w+ -J w .+pJ 2 z +2xJ z -J' + yJ' 2 ] +••• (18) 

V2 gfy sec 2 6 W 

where J z = J w ^—sm 2 6^J EM . The modified Weinberg angle, 9*, is introduced to account for 
the extra component of the electromagnetic current that mixing introduces. The experimental 
determination of the p parameter assumes that the non-oblique terms (x, y) are negligible. 
For the model discussed in Section [21 the non-oblique coefficients x, y presented below are 
determined at tree level, using eq ([2]): 

p sin ( sec x 



(cos C + sin 6 W tan x sin () 



(19) 



x 2 



y = — (20) 
P 

As | sin x | — > 1 (maximal kinetic mixing), |x| — > p esc W . In this limit, the product of the 
£7(1)' coupling and £7(1)' charges of SM matter must be very small in order to satisfy the 
experimental constraints. The tree-level p parameter in this model is: 
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Figure 1: Allowed parameter space for the Z' model given the measurement of the p parameter 
from electroweak precision data, with Q' H = 2/5 (the contours of p are spaced every 1CP 4 ). 
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(21) 



Once radiative effects in the SM are accounted for, the p (= p/p SM ) parameter, which is 
sensitive to physics beyond the SM, has been determined by a global fit to the EWPD giving 
p = 1.0002to;oooI with the la limits shown [25]. 

Fig Q] shows the allowed parameter space given this constraint (with Q' H = 2/5) for two 
values of a' = g' 2 /Att. From that one can easily extract the mass bounds on Z' once we know 
the values of the kinetic mixing, p parameter and a' . Given the definition of the p parameter 
in eq (JTHJ) , one sees that the limits presented are insensitive to the hidden valley at tree level 
and also to the SM fermion U(l)' charges. In Fig la, the point where the current coupling to 
the physical Z is identical to the SM case (C = 0) is at sinx ~ 0.25. The plot demonstrates 
that any value for m z , will satisfy the p constraint at this point when the tree level result is 
considered. For small m z ,, although not noticeable in Fig la, the 'arm' of allowed parameter 
space extending from the opposite sign of sinx does not meet with the straight p = 1 line 
(ie there is no 'hole' in the allowed parameter space). For large m z ,, there is an approximate 
symmetry in \ ~ > ~ X °f t ne allowed parameter space as suggested by eq (f2~T|) , 

As g' Q' H is increased, the position of the straight p = 1 line shifts to larger values for the 
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kinetic mixing. The 'arm' of the allowed parameter space also retracts, and the contours of 
p shift to larger values of m z ,. As g'Q' H is increased beyond sy/2, the picture of the allowed 
parameter space then resembles Fig lb and the contours continue to be shifted to larger 
values of m z ,. It should be noted that when the p constraint is very weak, other EWPD 
measurements provide stronger constraints. However, the strongest constraint on the model 
introduced in Section [2] is most often from either the p parameter or direct searches. 

3.2 Direct Z' searches 

In this section we consider the mass bounds produced by the Tevatron from searches for new 
gauge bosons beyond the SM. These can provide a more stringent lower bound for m z , than 
electroweak precision data in some cases, and is still present when there is no Z 1 — Z^ M mass 
mixing (£ = 0). Fig [2] shows their experimental limits [26], and the theoretical cross-section 
for pp —* Z' — > e + e~. 




Kinetic Mixing, Sin[;f] 



Figure 2: CDF limits on a new spin-1 Figure 3: CDF lower mass bound on the 

particle (FERMILAB-PUB-07-367-E) Z' in the model presented in Section EJ 

The boson labelled Z' has a lower mass bound of 822 GeV. Note that this Z' is the same as 

A, A, 

Z 1 introduced in Section [2.H but with no kinetic mixing and with the U{1)' coupling fixed 
[27] to g 1 = \J~^9y/ tan^. It is also assumed that the branching fraction to invisible states is 
negligible. Relaxing this assumption would reduce the theoretical cross section for producing 
SM final states, and so reduce this lower bound on the Z' mass. However, it will follow 
that this is a reasonable assumption for the hidden valley extension considered if the dark 
matter candidate discussed in Section is the lightest HV state coupling to the Z ', and the 
SM particles have similar (or greater) order of magnitude U(l)' charges compared to the HV 
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matter. The constraints on a Z' where the SM particles have no U{\)' charges have been 
discussed in other papers [28l l29| l30l 1511 [32] . 

The lower bound on the mass of the Z' introduced in Section 12.11 for arbitrary kinetic 
mixing and negligible decay to HV states is shown in Fig Since the theoretical resonant 
cross section is controlled by branching fractions, the lower bound on m z , only varies by 
~ 10% as the (g',x) parameters vary. For large a', the p constraint is much stronger. 

4 Dark matter constraints 

Let us now address the dark matter constraints on the model. We first state the method for 
calculating the thermal relic abundance of a stable particle and discuss what conditions are 
necessary for a dark matter candidate from the "hidden valley" to contribute significantly to 
the present day non-baryonic matter energy density. We present limits on the parameters of 
the hidden valley model from ensuring that dark matter is not overproduced and consider the 
constraints from direct dark matter searches. We find in particular upper mass bounds on Z' 
which will be combined with the lower mass bounds of Section [3j 

4.1 Thermal relic energy density 

During inflation, the reheating may have been sufficient for the Z' mediator to thermalise the 
SM and HV sectors. If not, the temperatures of the two sectors will be independent, and 
depend on the relative couplings to the inflaton. We consider the former case with a dark 
matter candidate from the hidden valley that annihilates into the SM sector. 

The thermal relic density of a stable particle is controlled by when this particle leaves 
thermal equilibrium with its annihilation products ("freeze-out"). This occurs roughly when 
the thermally averaged annihilation rate of this particle becomes slower than the expansion 
rate of the universe: 

n(cT atm v)~H (22) 

where n is the number density of the particle freezing out, (c ann v) is the thermally averaged 
annihilation cross section weighted by the relative velocity, and H is the Hubble constant. 
The ensemble of frozen-out particles then expand isenthropically, allowing the present number 
of particles in a comoving volume to be calculated. The current contribution to the energy 
density of the universe can then also be determined. If the particles are non-relativistic at 
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the time of freeze-out, the thermal average is given by: 



3/2 roc 

(a m v) « — = / dvv 2 (a ann v)e- xv / 4 (23) 



2vWo 

where x = m DM /T. The annihilation cross section can be calculated perturbatively, however 
for slow moving dark matter there can be important non-perturbative effects due to inter- 
actions before annihilation. This is accounted for by the Sommerfeld factor as discussed in 
detail in Appendix iDl and other papers [Ml ESJ [361 ETJ EH [39, 40J. The temperature at which 
the freeze-out occurs is determined from eq ([22]) . and is parametrised by the value of x at 
freeze-out [4*2] : 

x _ ln / [3.85 x 10 17 GeV] g DM m DM (fr ann v) Xp \ 

\ \/ x F 9*s J 

where g DM is the number of degrees of freedom of the dark matter candidate. The g* s pa- 
rameter is the effective number of relativistic degrees of freedom contributing to the entropy 
density of the universe. For T > 300 GeV, g* s = 106.75 (assuming only Standard Model 
particles contributing), and at T = 1 GeV, g* s ~ 80. The contribution to the energy density 
is then given by: 

_ , 2 8.6 x 10- 11 GeV~ 2 , , , f°° . _ 2/ . 

dm = / . \ — — where J(x F ) = dxx (cr ann u) (25) 

\Zg* s {x F )J{x F ) Jx F 

where the approximation that all annihilations cease after "freeze-out" has not been applied. 
The "annihilation integral," J(x F ), accounts for the reduction in particle number after "freeze- 
out." The thermal relic will contribute to the non-baryonic matter energy density observed in 
the universe today. The 5-year data from WMAP p] suggests that VL c h 2 = 0.1099 ± 0.0062, 
so the central value will be used as an upper limit for the dark matter contribution to Qh 2 . 

The thermal relic abundance increases as the annihilation cross section is reduced. Thus, 
whereas the constraints in Section E] led to small cross sections being preferred for Z' mediated 
or mixing events, the WMAP constraint will limit how small these cross sections are allowed 
to be, given a thermal relic from the hidden valley. 

4.2 Hidden valley dark matter candidate 

For the model discussed in Section "2j there are no flavour changing interactions for the HV 
fermions. A state with a non-zero "flavour" quantum number is therefore stable. If there are 
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confining interactions for the HV matter, some of the composite states will be flavour neutral 
(eg 7r°) and thus unstable due to the presence of the HV-SM mediator. However, there will 
also be stable states present (eg 7rJ) which are dark matter candidates. The lightest stable 
states dominate the relic density. Since the last annihilation mechanism to freeze-out controls 
the relic density, the freeze-out temperatures for the processes below must be evaluated and 
compared, to determine what the dominant annihilation mechanism is: 

• 7T+7T- -► Z' -» SM 

This process freezes-out at a temperature T x ~ m^± /25 when dark matter is not overpro- 
duced but the annihilation is not close to the Z' resonance (in which case 2\ ~ m w ± /40) . 
This was found using the annihilation cross sections to SM products listed in Ap- 
pendix O The analysis in the following section assumes that this is the dominant 
process at freeze-out. For 7r^ masses in the range 400 GeV to 40 TeV, the freeze-out 
temperature is passed at a time between 10~ 13 and 10 -8 sec. The Z' resonance is at 



4m^ ± ~ m Z i J m z , + T 2 Z , and this relates the upper mass bounds of the and Z' 
presented in Section 14.31 

7r+7r~ — > vr^7r° (followed by decay to SM products) 

This annihilation proceeds via the strong interactions of the hidden sector, and may 
freeze-out after the Z' mediated annihilation. Using the effective Lagrangian in eq (|12[) . 
the annihilation cross section is: 

36tt 3 Pf 

where the velocities of the initial and final states, j3 i f, account for the phase space. A 
Sommerfeld enhancement from v-strong or other interactions could significantly raise 
this cross section further. The freeze-out temperature can be found using eq 



771 ± 



a5 - ta (i^v)- 3ta fe) +ta W 



(27) 



The thermal average of the 7r^ velocity is found by integrating over the kinematically 
allowed phase space, which in the non-relativistic limit gives: 

Wf) - (V") e " M2)x/2 ^[( 1 -^ x/2 ] ' li6<1 (28) 
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where 5 = mjrrij, and K 1 is a modified Bessel function of the second kind. If the 
SM sector does leave thermal equilibrium with the HV before the v-strong interactions 
freeze-out, the total energy in the hidden sector is reduced by the decay of the flavour 
neutral v-pion to SM states j3J. The relic abundance of the 7r^ is then exponentially 
suppressed by the ratio of the time of freeze-out and 7r° lifetime, which can lead to a 
negligible contribution to the non-baryonic matter energy density. The lifetime is 
given by: 

• If the v-pions are nearly degenerate with or heavier than the Z' , the "valley" feature of 
the hidden sector is lost, but extra annihilation channels can become important for the 
7r^ dark matter candidate. This situation is beyond the scope of the paper. 

The unitarity limit places an upper bound on masses for a thermal relic at O(10 2 ) TeV 
At this limit, we see from eq (I27D that the v-strong interactions will freeze-out before 
the Z' mediated interactions if either A HV > 10 3 TeV or there is phase space suppression due 
to the 7Ty, 7r^" mass difference. If the 7r° is heavier than the 7r^, this process is known as a 
"forbidden" channel [45] . 

For TeV scale dark matter, we see that T 2 > T x if either A HV /m ± > 30 or there is phase 
space suppression, using eq ([2*7]) . In the case A HV /m ± = 8 (14), the Z' mediated annihilation 
will freeze-out last if the unstable v-hadrons are at least 10% (5%) heavier than the strongly 
interacting dark matter candidate. It has been shown in Section [2] that this kind of mass 
spectrum can easily be realised in the hidden sector and in the following we assume that 
T2 > T\ is indeed the case and the dark matter is 7r^. 

4.3 Allowed parameter space and combined mass bounds on Z' 

The constraints on a spin thermal relic which dominantly annihilates via the Z' are now 
considered. In this the annihilation cross sections of Appendix O are used, corrected by the 
Sommerfeld factors detailed in Appendix [*""". Figures [5] and [*"] show the allowed parameter 
space for the dark matter and Z' respectively, for certain values of f7(l)' coupling oi and 
charge Q, with the SM U(l)' charges as given in Table 1. These plots assume that either the 
dark matter mass or the Z' mass is a free parameter to scan over, with m DM < m z ,. If either 
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Figure 4: Allowed parameter space for the spin dark matter with U(l)' charge, Q, that 
dominantly annihilates via a heavier Z' with the SM U(l)' charges given in Section 

of these parameters are further constrained by experiment, the parameter space of the other 
quantity will also be further constrained. 

The upper mass bound on the Z' is found when the dark matter annihilation cross section 
on the Z' resonance is below a critical point, leading to over-production of the DM candidate 
with respect to the WMAP constraint. The lower mass bound of the dark matter allowed 
parameter space is controlled by the lower bound of the Z' mass and is the analogue of Lee- 
Weinberg bound for heavy neutrinos |46j . If the CDF constraint on the Z' mass is relaxed 
due to a significant branching fraction into HV states, it can be possible for the dark matter 
to be lighter than that shown in Fig 0] for small kinetic mixing (| sinxl ^ 0.4). 

As c/Qefi increases, where Q cS is the coupling of the to the Z' which accounts for kinetic 
and mass mixing using eq the effect of interactions before annihilation (the Sommerfeld 
effect) becomes much more significant. For spin dark matter, the repulsive interaction 
between the dark matter particles leads to an exponential suppression of the annihilation cross 
section. So, the allowed parameter space cannot be made arbitrarily large by increasing the 
U (1)' charge of the dark matter. When a' = 1 and the other parameters are unrestricted, the 
greatest value for m z , consistent with the WMAP determination of Q c h 2 is ~ 18 TeV. As the 
coupling of the dark matter to the physical Z' is proportional to secx, the Sommerfeld effect 
also produces much stronger limits on the amount of kinetic mixing allowed than expected 
from perturbative calculations. 

The upper limits on the masses of the spin dark matter and Z' are found where the 
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(a) a' = g' 2 /(4iT) = 1(T 3 (b) a' = </7(4vr) = 1 

Figure 5: Combined constraints on the Z' with the coloured regions showing the p parameter 
allowed space with la limits determined from electroweak precision data, and lines for the 
lower mass bound from Z' searches at CDF and upper mass bounds given the presence of a 
spin dark matter candidate with U(l)' charge, Q, that dominantly annihilates via a heavier 
Z'. (The contours of the p parameter are spaced every 10 -4 ). 

dark matter candidate saturates the non-baryonic matter energy density, J7 c /i 2 . If a separate 
particle is found at the LHC, which contributes to £l c h 2 , the allowed parameter space of the 
HV dark matter candidate and Z' will be reduced. However, if the annihilation cross section 
of the HV dark matter has been underestimated, the upper limits are raised. This could be 
the case from having neglected annihilation products involving supersymmetric SM states. 

It is also possible for the HV dark matter to dominantly annihilate through the physical 
Z, due to kinetic and mass mixing. This does not change the bounds on the physical Z' mass, 
but does sometimes allow a very narrow region of parameter space for the dark matter around 
the Z resonance, m DM ~ m z /2. This has not been shown in the plots. 

4.4 Direct dark matter searches 

This section considers the constraints from direct dark matter searches on the HV dark matter 
candidate. There are numerous experiments which have looked for recoils of nuclear matter 
due to collisions with dark matter. There have not been any signals observed, except for 
an annual modulation in the DAMA experiment which prefers light dark matter. The dark 
matter in this model is too heavy to explain the DAMA signal, so the upper limits on the 
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scattering cross sections from the other experiments are considered. The current and projected 
experimental limits are shown in Fig [6] with solid and dashed lines respectively. These cross 
section limits assume a local dark matter density of ~ 0.3 GeV cm -3 and mean velocity of 
~ 220 km s . If the dark matter candidate does not saturate the non-baryonic matter energy 
density, then the experimental limits are weakened. 

Fig [6] also shows the theoretical prediction from scans of the allowed parameter space, as 
given in Section f4.3l For the allowed dark matter masses, the Sommerfeld effect is negligible 
for these calculations. This is because the Z' is much heavier than a proton so the interaction 
is screened to a small distance compared to the characteristic scattering length scale. 

For a certain choice of parameters, the effective vector coupling of the proton to the Z' is 
zero. In the neighbourhood of this region, the dominant spin independent scattering occurs 
via the Z . Since the dark matter only gains milli-Z-charges from mixing, the cross section is 
suppressed. The vast majority of parameter space gives a theoretical cross section within a 
factor of 10 of the upper boundaries for the regions shaded in Fig El It will be unlikely that 
all experiments use atoms for the dark matter to scatter off that are in the special region of 
parameter space that allow the very small cross sections. 
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4.5 Indirect dark matter searches 

In this section, we briefly review the prospects for indirect detection of the HV dark matter 
candidate. There are searches for signals of dark matter annihilation or decay in our galactic 
neighbourhood. For the latter, the lifetime of the dark matter candidate should be greater 
than the lifetime of the universe in order for there to be a significant energy density currently 
present in the universe. Some possible signals include highly energetic photons, neutrinos 
and antimatter subject to a propagation dispersion relation. The stability of the HV dark 
matter candidate was inferred by the presence of a flavour symmetry. However, gravity is 
flavour blind and so we expect Planck suppressed operators to break this symmetry and 
mediate decay of the 7r^. As for dark matter annihilation signals, since the average velocity 
of cold dark matter in the galactic halo is of O(10 -3 ), the Sommerfeld suppression for the 
v-pions to annihilate via the Z' strongly reduces the cross section. This results in a negligible 
signal for current experimental detectors. For dark matter that aggregates inside stars, large 
planets or dwarf satellite galaxies where the velocities are O(10 -5 ), the annihilation signal is 
even further suppressed. Thus for the repulsive interaction considered here the prospects for 
indirect dark matter searches are poor, in contrast to the case considered in [39] that for an 
attractive interaction between the dark matter candidates, these objects can provide a much 
more favourable source for dark matter annihilation signals. 

5 Conclusions 

We considered experimental electroweak and dark matter constraints on a simple extension 
of the SM with an additional U(l)' massive gauge boson and a particular type of a hidden 
sector ("hidden valley") with a confining, QCD-like gauge group and hidden valley matter 
charged under the new U(l)' . As a result the SM and hidden valley sectors communicate 
via renormalisable operators involving Z' as a messenger, which can then probe the physics 
of the hidden valley sector. The latter can also provide a dark matter candidate. Combined 
electroweak and dark matter constraints placed both lower and upper bounds on the mass of 
Z' and these were studied in detail. 

We found from corrections to the p parameter in electroweak precision data and direct 
searches a lower limit on the Z' mediator mass of 0(1-10 TeV) depending on the U (1)' coupling 
constant and kinetic mixing parameter. This result followed when there is a negligible Z' 
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partial decay rate to HV states and 0(1) U{\)' charges for the SM matter. This limit does 
not require the presence of a hidden valley, as it would similarly apply for a pure U(l)' 
extension to the SM. In addition, it was demonstrated that the kinetic mixing of the U(l)' 
and hypercharge gauge bosons is not restricted to be small by current observations. 

We also found a dark matter candidate in the strongly interacting hidden valley that could 
saturate the present day non-baryonic matter energy density. If the spin HV thermal relic 
dominantly annihilates via the Z', and the U(l)' charges of the SM and HV states are of 
similar magnitude, there is an upper bound on the Z' mass of 0(10) TeV and restrictions 
on the amount of kinetic mixing allowed. Although the upper limit is beyond the reach of 
the LHC [29} [33] . a large proportion of the parameter space will be tested. Detection of an 
alternative candidate to contribute to the dark matter density at the LHC could forbid the 
presence of this HV thermal relic. The lightest HV state was found to annihilate close to the 
Z' resonance, restricting it to have an O(TeV) mass. Future direct dark matter searches will 
also significantly probe the currently allowed dark matter parameter space. 

The inclusion of the Sommerfeld effect for p-wave annihilations of the spin dark matter 
via a spin 1 state brought a significant change to the annihilation cross section. This con- 
strained the parameter space much more than perturbative calculations would have suggested. 

If the lightest HV states are O(TeV) and the model discussed is realised, then extra signals 
containing highly energetic SM decay products of HV states may be observed in collider 
detectors. The pseudoscalar 7r° state would preferentially decay to top quarks, gauge bosons 
and Higgs fields, and the LHC / Tevatron would see low multiplicity events. There would 
be insufficient energy and resolution to see the influence of the non-pseudo-Goldstone boson 
states. In alternative hidden valley scenarios, the phenomenology may be different and further 
study is important to identify how such models and their collider signals can be constrained. 
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6 Appendix 



A GUT embedding and neutrino masses 

We provide here further details to the discussion in Section [2] to show how neutrino masses can 
be generated in the model considered. The model discussed may be considered a low energy 
realisation of a theory with SO(W) symmetry if the SM Higgs field is part of a 10, 120 or 
126 representation of SO(W), and the non-SM states have quantised U(l)' charges, q±, Q'^ = 
0, ±1, or ±2, also forming parts of SO(10) multiplets. The missing states of these SO(10) 
multiplets can have masses of the 50(10) unification scale depending on the mechanism that 
breaks this symmetry. 

To generate a large Major ana mass for the right-handed neutrino via the Higgs mechanism 
we choose Q'^ = 2, and so if unification is desired <j> must belong to a 126 representation, 
which may be the multiplet H is embedded in. The vacuum expectation value of (f) is also 
responsible for generating a mass for the Z' which we will find is O(TeV) scale. The neutrino 
mass matrix for one SM family is of the form: 



The relation of the scalar vevs to gauge boson masses and couplings is given in Appendix O 
For A of O(l), the heavier neutrino has a mass similar to or heavier than the Z' boson, and for 
a phenomenologically acceptable mass of 0(eV) for the light neutrino, the Yukawa coupling 
leading to the Dirac mass term must be small, y u ~ 10~ 5 . This is not much smaller than 
the electron Yukawa coupling. However, if one wants to avoid the introduction of such a 
coupling, a Z 2 symmetry may be introduced as suggested in [3] with only the N i fields having 
odd parity (in this case the Ni fields are dark matter candidates). This assignment is discrete 
anomaly free if heavier states with odd parity have been integrated out [10], and it forbids 
the Dirac mass term. A Majorana mass could be generated for the left-handed neutrinos by 
a G SM x U(V) invariant dimension 6 operator, LLHHcf)* /M+ . If this operator is responsible 
for the light neutrino masses, it suggests that the scale of new physics which leads to lepton 
number violation is ~ 10 8 GeV. Figure lA-il shows some mechanisms for generating this 
dimension 6 operator. 




(A-1) 
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Figure A-i: Some Feynman diagrams that lead to an effective Majorana mass for the SM 
left-handed neutrinos after integrating out heavy scalar states, The "X,Y" label on & X Y 
indicates that this scalar can be considered as part of an X or Y representation of SO(10). 
The $ a field has the representation (1,1)(1,-|) under SU(3) x SU(2) x U(l) y x U(l) x , $ 6 
is (1, 1)(-1, -§), and $ c is (1, 2)(-i, f ). 

B Gauge boson eigenstates 

Here we provide additional details to the discussion in Section 12.11 by computing the gauge 
boson eigenstate fields in the presence of various mixings present. The kinetic mixing in 

£gau gc kinetic = ~J [ F V F Y + F ' ' F ' + 2 sin X F '' F Y ) + non-Abelian kinetic terms. (B-l) 
is removed by the redefinition: 




)l 






\ % ) 



(B-2) 



where is the SM hypercharge boson and B' is the U{\)' gauge boson. After this, B^ and 
B' have canonical kinetic terms. Under this transformation (b = C b) the source terms for 
the gauge fields in the Lagrangian, h T j, must be invariant, and this fixes the relation between 
new ("tilded") and original currents: j = C T j. It then follows that with the field redefinition 
in eq (|B-2p . the couplings to B^ (the new hypercharge) are the same as in the original basis. 
After EWSB it is convenient to introduce the following fields: 
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The field will have no mass terms since it couples to the generator T 3 + Y by definition of 
the Weinberg angle (U(1)em unbroken). The SM Higgs field couples to both the SM neutral 
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gauge boson and to B' , so after spontaneous symmetry breaking there will be mass mixing 
terms. In general, mass mixing terms still exist after going to the "tilde" basis: 

' % Z') I 7 I I ~. I +^ + ^-term. (B-4) 



/ m 2 -A i 




f z» \ 


I -A M 2 j 


1 





2 v M M 

The rotation necessary to diagonalise this mass matrix is parametrised by £ as follows: 
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(B-5) 



where Z± are the mass and kinetic eigenstates, called the physical Z and Z' in the main text. 
With these relations one finds eq.([2]) in the text for the transformation of the currents. 

C Cross section formulae 

We list here the partial decay rates and cross section formulae used to generate the limits 
presented in the text, Section [H The interaction terms considered are parametrised below, 
with canonically normalised kinetic terms assumed: 

iQxQx <P W) X" + 9 XX , ^ X»Xl 4> - 9x^1^ (c v - c a75 ) X 1> + kc. (C-l) 



where g xx , = g x g x , Q x {<p) Q x , {<f>) (C-2) 

where ^ is a spin field, = V2(<p}, X^,X'^ are gauge fields and g x ,g x , are the respective 
coupling constants. Q x is the charge of the (ft field under the symmetry whose gauge boson 
is X^ and so on. The decay width of the Z' is an important quantity for any Z' mediated 
event close to the propagator pole. Some partial widths are given below, with the others 
assumed to be negligible. When relevant, the final state spin or polarisation configurations 
are summed over (but not for different colours). 
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where W 3 ^ = c x X» + • • • ; (V 3 ^ = sin 9 W A» + cos ( cos 6 W + (- sin C cos 6 W ) z£) 
C.l Spin dark matter 

The parameters associated with the initial bosons are labelled (i), and those with the final 
states labelled (f). The tree level cross sections for annihilation via a vector boson are given 
to leading order in and are all p-wave annihilations: 

• Dirac fermion anti-fermion product 
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• Spin 0, Spin product 
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Gauge boson (X') and a scalar (<f>) product 
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D Sommerfeld effect 

In the following we provide more details about the Sommerfeld effect mentioned in Section |H 
A more complete version of this study will appear elsewhere [47j . When the initial or final 
state of an interaction event involves slow moving particles, attractive or repulsive forces 
between these particles can lead to large enhancements or suppressions in the cross section 
compared to a perturbatively calculated result. This has long been known in nuclear physics 
and condensed matter physics as the Sommerfeld effect. Recently, its importance has been 
pointed out for cosmological calculations of thermal relics, and the indirect detection of dark 

matter [HI E51 EH E3 EHl ESI HQ] . 

The non-perturbative physics can be thought of as the limit of perturbative Feynman 
diagrams, with an infinite number of particle exchanges. For 2-body scattering this involves 
a "ladder" diagram. The solution of the non-perturbative vertex function, T, is given by the 
Bethe-Salpeter equation: 

Pi Ps , 

P. \ ^ \ / pi+p 2 -<i 

Figure D-i: Bethe-Salpeter equation in diagrammatic form. 
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/d 
-j—^fG{q)G{p 1 +p 2 -q)T (D-l) 

where T are "compact" vertices which do not involve any intermediate state composed solely 
of the scattering particles, and G is the non-perturbative propagator. The inhomogeneous 
integral equation for T can be solved to arbitrary accuracy, given that T and the propagator 
are sufficiently well known from perturbative calculations. The finite ladders can be neglected 
when near the poles of the scattering states, which leads to the non-perturbative nature of 
this effect and a homogeneous integral equation. For an annihilation process as shown in 
Fig ID-iil the matrix element for the diagram with the ladder is related to the one without a 
ladder as given below: 

Z Final States 




Figure D-ii: Diagram of Non-Perturbative Scattering before Annihilation 



M with ladder (p 1 ,P 2 ;{Pf}) = J ^-j4 ■A^ w/o (p 3 ,P4;{p / })G(p3)G(p4)r(p 1 ,p 2 ;p3,p 4 ) 

^3 ^ w/o (k,-k;{ P/ }) V> BS (k), (D-2) 

where the Bethe-Salpeter wavefunction, ^ BS , is introduced in the centre of momentum frame, 
and the intermediate scattering particles are taken on-shell. It follows that for non-relativistic 
processes, the Bethe-Salpeter wavefunction is equivalent to the solution of the Schrodinger 
equation with a potential that accounts for the scattering interactions [H]. For central po- 
tentials, the spherical harmonic basis is convenient for solving the Schrodinger equation: 

^Elm( r ) = REl{ r ) Y lm( T,4>r) ( D "3) 

The Fourier transform of eq ()D-3|) is determined below: 
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It is therefore useful to decompose the matrix element in the same orthonormal momentum 
space basis to calculate eq (|D-2p . The weighting of the components of the Bethe-Salpeter 
wavefunction with different (l,m) quantum numbers is determined by the matrix element. If 
.M w/o is expressed as a polynomial in k, the result of integration of the wavefunction weighted 
by extra powers of k needs to be found. This can be done by considering derivatives of the 
position space radial wavefunction and the inverse Fourier transform relation: 
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Using a series expansion of the spherical Bessel function, then 
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Suppose the full matrix element is dominated by the contribution of one term in M v/o oc 
k l+2n Y lm , where n is some non- negative integer. In this case, the cross section which accounts 
for the ladder is related by an overall factor to the cross section calculated from diagrams 
without the ladder (<7 with laddcr = S ln c-), valid in the non-relativistic limit. S l n is called the 
Sommerfeld factor, with the value: 
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It is common to ignore the effect present for non-zero angular momentum processes, arguing 
that the angular momentum barrier always suppresses higher partial wave processes. However, 
if s-wave interactions are small compared to higher partial wave processes, as when scalars 
annihilate via a heavy vector, this argument is not always applicable. 
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D.l Coulomb potential 

The reduced Schrodinger radial equation for a two-body system of particles each with mass 
M, and a Coulomb potential is as follows: 



(D-9) 



where (3 is the speed of each particle when at infinite separation in the centre of mass frame. 
This can be re-written in terms of dimensionless variables: 

dz + T + 7T- ~ ^4-^ ) zR,(z;x) = where z = 2rM(5 x = AI(3 (D-10) 
4 2 z z A J 

The solutions are given in terms of a confluent hypergeometric function: 
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The Sommerfeld factor written in this dimensionless formulation is: 
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The result for non-zero angular momentum states due to Coulomb interactions is thus: 
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In a perturbative expansion, higher partial waves are suppressed by factors of j3 2 . However, 
in the limit of small velocities (large x) the Sommerfeld factor for a Coulomb interaction fixes 
each partial wave to have the same velocity dependence. 



D.2 Yukawa potential 

The Schrodinger equation with a Yukawa potential is related to the Coulomb case as follows: 



Ae 



Ae -xz/(2y) where y = AM / m 
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(D-15) 



where m x is the mass of the mediator. The Sommerfeld factor for the / = case has been 
presented by Cirelli, Strumia and Tamburini [38]. Figure D-iii shows the I = 1 case (n=0), 
where numerical simulations have been used to find the wavefunction. 
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Figure D-iii: Sommerfeld Factor for Yukawa interaction of a p-wave state (contours not 
labelled vary by a factor of 10) 



For small y, the Sommerfeld factor gets closer to unity as the angular momentum increases. 
This is because as the interaction is screened beyond some short distance, the increased angular 
momentum barrier becomes more efficient at keeping the wavefunction away from this core. 
For a repulsive Yukawa interaction, there is then less suppression for the higher partial waves. 
So, although higher partial waves may be neglected at a perturbative level, this ladder effect 
could cause them to become significant or even dominant. However, for attractive interactions, 
the larger angular momentum processes are enhanced less, so higher partial wave terms which 
are negligible in the perturbative expansion stay negligible. 

There is always some reduction in the velocity dependent factor of the partial wave cross 
sections, in going to higher partial waves for non-zero (3 , but the P 2Al suppression is only 
found for large velocities. For the attractive Yukawa potential, / = 1 bound states exist when 
y > 9.08. The system can be close to the Breit-Wigner tails of these resonances when the 
relative velocity is small, and the Sommerfeld factor reflects this. In the limit y — > oo, the 
Coulomb potential is approached, so the discussion in the previous subsection is relevant. 
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